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Hn ,T-TxmAgTVP r. ATmTAC OUT Pttt AMD PULMONARY FUNCTION 
MONITORING T?RSPTR ^T> OAR ANALYSTSTECHNIQU^ 
AND PHYSIO 1 nmr. AT MODELING 

PKOSS REFEF™fF TO RET ATFX) APPLICATIONS 
This application claims priority trom U.S. Provisional Patent Application Serial No. 
60/116,648 entitled "Non-Invasive Cardiac Output Monitor Using Respired Gas Analysis 
Techniques And Physiological Modeling", filed January 21, 1999, and from U.S. Provisional 
Patent Application Serial No. 60/140,763 entitled "Non-Invasive Cardiac Output and Pulmonary 
Funcuon Monitoring Using Respired Gas Analysis Techniques And Physiological Modeling* , 
filed June24, 1999. The disclosures of these proviaonalapph^ 
reference in their entireties. 



R ACKGROUN n OF THE INV 



fi^A nf the Tnvention 

The present invention relates to methods and apparatus for determining cardiac output, 
,he amount of blood the heart is pumping, as well as identifying pulmonary functions, without 
resorting to invasive techniques which introduce foreign objects and me like, into a body. 



TWrri ption o f the Related Art 

Monitoring me cardiovascular system to determine myocardial performance is of 
paramount importance in patient care, regardless of whether the patient is located in the 
20 physician's office, emergency OT operating room, inte^^ 

transit (eg., in an ambulance). Although routine cardiac monitoring usually begins with a 
determination of the patient's heart/pulse rate and blood pressure, in the case of patients who are 
eq^encing cardiac difficulties or distress, additional diagnostic details regarding the operation 
of the heart are needed. Such additional monitoring may quickly progress to include an 
electrocardiogram (ERG) and the measurement of hemodynamic variables such as cardiac output 
The term cardiac output is defined as the mean or average total blood flow in the 
ciculaiorysvstmper^^^ Cardiac outout is assodatrf^ 

is consequently an important parameter in assessing the condition of a patients health. 
Knowledge of cardiac output level and trends have important diagnostic value in that mey 
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provide the clinician with information to help him/her assess how well me myocardium is 
functioning so as to provide the basis for the timely delivery and prescription of appropriate 
therapeutic modalities. Pulmonary function relates to the ability of the -body to make use of 
oxygen and to eliminate wastes such as carbon dioxide, this parameter is strongly affectedby 
5 physiological conditions such as deadspace and shunts, and the ability to quantify these 
conditions forms a major part of cardiopulmonary therapy. 

Owing to the uncertainties of the geometry ofblocd vessels (e^ diameter, comply 
etc.) and the dynamic nature of the heart itself, conventional flowmetering techniques such as 
flow resistance measurement or velocity (e.&, Doppler and ultrasonic) measurements have 

10 proven unreliable in estimating cardiac output As a result, cardiac output is routinely measured 
invasively; that is, by surgically placing an instrument into the arteries near the heart 

The current state-of-the-art, and arguably the -gold standard" for cardiac output 
measurement, is considered to be either the Direct Fick or mermodOution technique using a flow- 
directed catheter (Swan-Ganz catheter). The catheter is physically threaded through a large vein 

1 5 (femoral, internal jugular, etc.) into and through the right atrium and right ventricle of the heart 
into fcepulmc^ artery 1^ At tbirt point, Ihem^ dilution 

techniques may be used to quantify the blood flow. Urifortunately, because of me mvasrvenature 
of the technique, the potential risk to the patient of hemorrhage, dysrhyumua or cardiac arrest is 
relativeh/high. Consequently, the routine use of mvasivetechmquessuchasmenn^ 

20 measure cardiac output is presentiy h^ted to specific clinic^ 
outweigh the risks. 

A significant number of patients (as many as two percent) do not survive me surgery 
associated wim^ 

situations where patients are extremely ffl and me m 

25 mortality is acceptable. Efficacy studies in recent medical literature report data mat raises 
questions as to the risk-benefit ratio of the information provided by invasive cardiac output 
measurement and whether invasive cardiac cutout measurement is in the best inteiest of the 
patient In addition to the mtrinsic danger of the invasive procedure, the monetary cost of the 
procedure is relatively high, as it is m itse^^ 

30 procedures,itshanaVonlabormtensityb^ ft 
is estimated that nearly $200 million is spent mtte U.S. alone for mvarivecaro^ 
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equipment and materials, and recent medical literature also questions the cost effectiveness of 

these invasive techniques. 

Adding to the increasing concern that the cost-benefit ratios may not be in the patient's 
best interest is the fact that many patients may not be adequately monitored and consequently are 

5 being put at risk from lack of diagnostic information. Presently, no reliable, accepted, cost* 
effective, non-invasive techniques are available for continuous monitoring; thus, only the sickest 
and highest risk patients are candidates for continuous cardiac monitoring. This leaves a huge 
population that goes unmonitored, of which it is well known that significant numbers encounter 
cardiac distressed Therearemany 

10 ciirical situations 

and home care where monitoring cardiac output is not routine, but if it were, would be of 
sigrJficar^benefittopatientcare. Therearesign^ 

n^any of which may have been prevented had myocardial function monitoring been available and 
appropriate responses initiated. Current estirruues o^ ^ 
15 cardiac compUcations exceed $22 bulionin the U.S.alonc 

Non-invasive and less invasive techniques are merefbre highly de^ 
because of the variability and complexity of me physiology of me circulatory svsteman^ 
pathology of disease, no c^ 

to be capable of obtaining reliable cardiac output values. Although less-invasive methodology 
20 such as impedance cardiography, Dc^pler-s^ 
singlo-br^ 

i„ their current implementations, they all suffer from significant problems and/or disadvantages. 
In general aU of mese techmques are ext^ 

staff, and are limited to a few well-defined clinical situations. In addition, each technique has 

25 unique specific lrmitations. ^ , V . 

More particularly, impedance cardiography requires the correct placement of electrodes 

on the neck and abdomen that are excited by a high frequency (e^ 00^ 
subsequent monitoring of the resulting impedance changes between the electrodes. The 
tapedance changes of the chest are used to determine the cardiac stroke volume resultmg from 
30 the expansion and contraction of the cardiac volume. Cardiac output can be calculated by 
combining this volume with heart rate in an appropriate algorithm. The limitations of this 
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technique include: the need/ability to correctly place the electrodes, accurate accounting for the 
volume changes resulting from the inhalation and exhalation of the lungs, and patient movement 
Furthermore, the high impedance electrodes act as antennas that pick up considerable amounts 
of electromagnetic interference (EMI), thereby interfering with the measurements. 
5 The Doppler-shift technique is based on the effect of the shift in frequency of sound from 

a stationary source that is reflected by a moving object With this method, the average velocity 
of the blood flowing in an artery can be readily measured However, to determine the volumetric 
flowrate, the cross-sectional area of the artery must be known. Obviously, soft tissue 
visualization techniques such as MRI are not practical at this time for general use, and ultrasound 

10 imaging generally tends not to be accurate enough, although it is used to provide a relative 
measure in some applications such as transesophageal-echocardiography. Costs are prohibitively 
high, and in this age of managed care cannot be considered practical for routine use. Esophageal 
Doppler techniques are also plagued with inevitable patient motion artifacts. 

A variety of indirect Fick techniques, including breath holding, single breath and 

15 rebreathing, have been proposed over the years to estimate cardiac output from various 
measurements of respired and tracer gasses. Breath holding and single-breath techniques using 
tracer gasses have had limited success but are not suited to continuous monitoring. A single- 
breath technique proposed by Kim et al. in "Estimation of true venous and arterial PCOz By Gas 
Analysis of a Single Breath,- J. AppL PhysioL, Vol. 21, No. 4, pp. 1338-1344, (1966), 

20 incorporated herein by reference in its entirety, probably had the greatest potential because of the 
promise of breath-by-breath monitoring. However, general acceptance has been lacking due to 
technology limitations for precise, real time, simultaneous respiratory gas measurements and 
limited experimental validation of their underlying assumptions. 

Perhaps most popular and most widely accepted of the indirect Fick techniques have been 

25 CO2 rebreathing techniques. Presently, the only commercially available noninvasive device is 
the Novametrix Non-Invasive Cardiac Output (NICO) monitor, which monitors respired carbon 
dioxide production combined with partial rebreathing (inhaling air with elevated carbon dioxide 
levels) in a variation on the well-known Fick Principle. (The Fick Principle relates essentially 
to a statement of flow continuity and mass balance over the cardiovascular system.) Mote 

30 specfficaMy, a non-dispersive infrared 

concentration and respired volumetric flowrate and hence CO2 production. The Fick equation 
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is used to calculate cardiac output as the ratio of the carbon dioxide produced to the arteriovenous 
difference of carbon dioxide content inblood. NICO is reported to have reasonable correlation 
with direct Fick and indicator dilution measurements in patients with normal, healthy lungs with 
minimal deadspace and/or no pulmonary shunts. 

However, the NICO system's reliance on the products of metabolism (U., the 
Novametrix sensors can measure only carbon dioxide) results in questionable accuracy in the 
presence of shunts and deadspace in the lungs; the accuracy of the NICO system is also 
compromised because it must rely on compensatory algorithms that are highly dependent on 
physiological conditions and unknown metabolic and respiratory parameters. Consequently, 
results are poor for patients with pulmonary and/or obstructive airway disease due to the effects 
of V/Q mismatching caused by increased pulmonary shunts and deadspace. These effects 
invalidate the assumption that PeCOj can be used to approximate the values for PyCOj and 
P.CO* The dilutional effects of a significant shunt on the pulmonary capillary blood flow 
invalidate the assumption that systemic cardiac output is equal to pulmonary capillary blood 
15 flow. Reasonable success has been achieved in compensating tor shunts by measuri^ 

of Ch saturation in a peripheral artery with a pulse oximeter. The major disadvantages of mis 
technique are: bulk of the rebreathmg ar^aratus and the time re^ 

calculate cardiac output This latter disadvantage precludes the use of this device for craitinuous, 
or even breath-by-breath monitoring; consequently, dynamic changes may not be detected quickly 

20 enoughforprevenuvemeasura^ 

man a recirculation time, readmgs may be affected by accumulated COa. 

While monitors for the continuous, breath-by-breath, measurement of CO* O* and 
anesthetic agents are commercially available, all are lacking in one or all of me following 
attributes: rebability, ease of operation, accuracy, the need for calibration, small size, and low 

25 acquisition cost and life-cycle cost For example, CCfe monitors using non-dispersive m. 
spectroscopy can cost over $1000 for hand-held versions and as much as $20,000 (with 
additional high life-cycle costs associated with the periodic calibration and maintenance of me 
equipment) for a full-spectrum operating room gas monitoring. Mass spectroscopy and Raman 
scattering systems are even more costly and bulky. In addition to cost, physical size, and 

30 inconvenience of operation (cahbration) conventional systems have found limited use of gas 
monitoring in the field for such things as validation of endotracheal (FT) tube placement during 
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emergency intubation and patient transport. Extubation, leading to severe, irreversible 
consequences, frequently occurs during patient transport, yet no monitors meeting the above 
characteristics have been available. 

Consequently, there remains a need for a reliable, cost-effective, non-mvadve cardiac 
5 output monitoring system capable of continuously measuring cardiac output and pulmonary 
function on a breath-by-breath basis usmg measurements of insr^ The 
availability of cardiac output measurement to routinely monitor the large population currently 
without benefit of such monitoring could significantly reduce the huge aftercare costs and 
morbidity and mortality resulting from undiagnosed cardiac complications in rion-cardiac-related 

10 procedures. A Ughtweight, rugged device would be ideally suited for use in field environments 
such as the ambulance and MEDEVAC transport, as well as the doctor's office, clinic, emergency 
and operating rooms and in intensive care units (ICU). 

Since there is no currently acceptable noninvasive cardiac output monitor available for 
routine use, mere remains a need for a technique to accuiatdy measure cardie 

15 eliminate risk of infection or invasive trauma to the patient Further, any technique that is 
economical, reliable, accurate, and simple to operate and maintain becomes a candidate for 
routine utilization. Moreover, a device that is Ughtweight and smau, opem me market to 
ambulatory monitoring, sports and physical fitness, and home care of cardiac patients. Finally, 
such a device would complement rural and military telemedicine where remotely located 

20 specialists can diagnose and treat patients given sufficient patient data. 

STTMMARV OF JW INVENTION 
Therefore, m Kght of me above, and for omer reasom 
invention is fully described, an object of the present invention is to provide a non-invasive 
cardiac cutout monitoring system that uses measurements of inspired and respir^ 

25 determination of cardiac output 

It is anomer object of the present invention to utilize a mathematical model of me human 
physiology that will compensate for variations in physical and disease states in determining 
cardiac output from respired gasses. 



WO 00/42908 PCT/USOO/01465 

7 

It is a further object of the present invention to measure uptake and release of inert and/or 
insoluble indicator gasses that are not metabolized and absorbed in order to eliminate the vagaries 
of the metabolic and absorption processes in determining cardiac outjput 

It is yet a further object of the present invention to use a gas analyzer that measures or 
5 assays all the gasses that are inhaled and respired, not just an indicator gas alone, thereby 
allowing for a complete description of the uptake, distribution and release of the gasses, mat then 
allows for accurate inputs to the physiological model. 

It is still a further object of the present invention to provide for a very low cost 
implementation of the technology in order to promote widespread use and to improve the general 

10 standards of care for patients. 

It is another object of the present invention to measure, in a real time, breath-by-breath 
situation, oxygen and carbon dioxide concentration from which both mixed venous and arterial 
concentrations of carbon dioxide can be determined. 

Another object of the present invention is to measure on a real time, breath-by-breath 
15 basis the anatomical and physiological deadspace of the lungs by combining breaming mass flow 
measurement with concentration waveform analysis. 

Still another object of me present invention is to provide a low cost means for 
determining the cardiac output and pulmonary function of a human bemg on a breath 
basis while accurately accounting for disease states as well as physical conditions. 

Yet another object of the present invention is to provide a cardiac output monitoring 
device that measures attributes of respired gasses on a bream-by-bream basis, which 
measurements can be used with any of the krtownFick techniques for determining cardiac output 
non-invasivcty. 

A fundamental aspect of the present invention is the use of a respired gas analyzer mat 
25 is capable of simultaneously quantifying the concentrations of several gasses in real time and a 
true, real time mass flowmeter to calculate uptake, production and expiration of gasses to provide 
measurements with known relationships to cardiac output and pulmonary function. A gas 
analyzer suitable for use in the present invention is disclosed in pending U.S. Patent Application 
Serial No, 09/104,997 entitled "Method and Apparatus For Real Time Gas Analysis" filed June 
30 26, 1998 by Tadeusz M Drzewiecki, and a pending provisional U.S. Patent Application Serial 
No. 60/121,370 entitled "Methods and Apparatus for Real Time Fluid Analysis'' filed February 



20 



wooo/tt«» rctnantmms 

8 .. 1 

25, 1999 by the same inventor. The subject matter disclosed in those applications is incorporated 

herein by reference in its entirety. 

The combination of a real time mass flowmeter and an inexpensive gas analyzer capable 
of simultaneously determining concentrations of multiple gasses in real time permits for the first 
5 time accurate determination of cardiac output on a bream-by-breath basis from analysis of 
respired gasses. More particularly, the cardiac output monitoring system of the present invention 
can be used with any of the Fick-principle-based non-invasive techniques that have been 
proposed in the art for measuring cardiac output from respired gasses, but that have heretofore 
been impractical, prohibitively expensive, inaccurate and/or unreliable. Moreover, the 

10 parameters measured and the extensive information provided in real time by the cardiac output 
monitoring system of the present invention allow known techniques to be refined and extended 
to more accurately account for pulmonary factors such as shunts and deadspace in the 
determination of cardiac output 

The gas analyzer disclosed in me aforementioned Drzewiecki patent applications 

15 simultaneously and in real time assays gasses, allowing accurate quantification of all the 
constituents of respiratory gas mixtures. Because me gas anaryzer measures 
of a gas mixture, including density and viscosity, a conventional flowmeter can be compensated 
for changes in gas properties, not only as a function of temperature but also for changes in 
composition. This allows the use of any one of a variety of low cost pressure-drop-type (fixed 

20 or variable orifice) flowmeters to accurately measure respired flows over a wide range of gas 
compositions, with equivalent accuracy of expensive mass flowmeters. Thus, artifacts caused 
by breathing in products of combustion or other gasses (e.g., anesthetics) that would affect the 
computation of gas uptake/production by giving erroneous volumetric or mass flows are 
eliminated. In this manner, me volumetric gain or loss from the lungs can be quantified 

25 throughout the respiratory cycle to provide me necessary d^ to accurately calculi 

output using the Fick Principle. By overcoming the technical difficulty of measuring the 
concentrations of oxygen and carbon dioxide simultaneously (with standard errors mat cancel 
rather than add as they do with independent sensors) an accurate measure of cardiac output can 
be obtained using single-bream techniques, such as that disclosed by Kim. Furraenndre, by 

30 providing an improved methodology for estimating alveolar COa and O2 concentration values 
that includes the effects of physiological (including alveolar) deadspace, by using the Bohr 
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equation combined with the considerable work of Fletcher on analyzing CXVvolume waveforms, 
in combination with an iterative anatomical/physiological model of gas exchange that converges 
on measured expiratory gas concentrations, the Kim technique provides accurate results under 
significantly broader conditions to include exercise and disease states. Finally, by including a 
5 pulse oximeter to measure Ch saturation, pulmonary shunts can be compensated for directly ia 
the expression derived for the Ch tension rather than by trying to estimate a value for shunts 
directly and thereby adjust the CCh values. 

The cardiac output monitor of the present invention poses essentially no risk to the 
patient, is easy to use, is inexpensive to manufacture and has virtually no low life-cycle costs 
10 (eg., no recaKbration is ever required), thereby making it economical to operate, andcanbe sized 
and packaged to be handheld while ^ mamtainmg an instrument (e.g., wave 
capability. 

By measuring physical properties such as density, viscosity and specific heat with very 
simple but highly precise pressure, flow, temperature and frequency transducers, the assay of the 

15 constituent concentrations is precisely calculated. The unique combination of concentrations that 
make up a gas niixture with given measured properties (visa>rity, denshy, specific heat) is 
determined by deconvolving the fundamental relationships mat define mixture property values 
in terms of their constituent concentrations. State-of-the-art, low cost, ultra-high dynamic range, 
microelectromechanical system (MEMS) pressure transducers, and a highly precise platinum 

20 RTD temperature sensor integrated with a specially designed fluidic oscillator flowmeter, 
measure the pressure drop, temperature and flow in a microfluidic capillary viscometer, orifice 
densitometer, and sonic microcalorimeter (specific heat sensor) integrated on a precision nucro- 
injection molded I^ratory^n-a-Chip (LOAC). A high-speed microprocessor provides 
solutions to the governing equations and drives an LCD. Because the concentrations ate 

25 determined from physical first prmdples, me gas analyzes ^ 
maintenance, which is a major advantage for field-use devices. 

According to one embodiment of me preset a known amount of an essentially 

inert, insoluble, indicator gas is tracked during respiration. The mput parameters (cardiac output, 
deadspace, shunts) to a validated physiological model (software) are iterated to obtain a matching 

30 time-history of the released (exhaled) indicator gas. The use of a model that allows for individual 
variations and disease states, as well as variations in body mass and uptake and distribution 
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parameters, results in a credible as well as accurate output, that, because it matches the measured 
values, represents a reasor^le estimate of the parame^ Hie values of cardiac output, shunts 
and deadspace that match the measured values are thus among the outputs of the monitoring 
system of the present invention. 
5 In another, more general embodiment of the present invention, the use of a tracer gas is 

dispensed with and the analysis of the consumed oxygen and produced carbon dioxide and their 
concentration waveforms is used to provide a measure of the mixed venous and arterial 
concentrations of carbon dioxide ar«U measure of me anat^^ 

leaving the physiological model to be used only to correct for and measure pulmonary shunts 

10 corresponding to a particular disease state. 

The low cost, affordable, accurate respired gas arialvsis technology that constit^ 
basis of the present invention provides a mechanism for determining the concentrations of the 
constituents of a gas mixture by measurement of certain independent physical and/or 
m ennodynarnic properties such as density, viscosity, specificheat, dielectric constant, refractive 

15 index, electromagnetic radiation absorptivity, etc., of the mixture and determming the assay of 
themixturethattMoa^ 

applied to the measurement of cardiac output, offers rigrimcain cost and diagnosis 

over other technologies that utilize the well-known and accepted Fick principle. For example, 

currently available non-invasive memods (e.g., me aforementioned 

20 the use of only one indicator gas (carbon dioxide or oxygen) at a time. Since these gasses are 
products of metabolism or are themselves metabolized, the algorithms used are necessarily 
complex andnotweuvahdated because mey must be able to accurately conritoti^ 
process. By being able to use essentially inert (non-metabolbdng) gasses, nitrogen and/or 
anesthetic gasses may be used as indicators and, moreover, me presence of more tiianone^ 

25 be monitored simultaneously. This enables the clinician to select the indicators), or 
combinations thereof, appropriate to a particular clinical situation. Forexar^le.adenrtrogenated 
patient on pure oxygen (and anesthetic gasses) in the operating room (OR) is a good candidate 
for a nitrogen indicator, whereas, a patient in the intensive care unit (ICU), who is breathing air 
or air and oxygen, may be a candidate for an anesthetic agent or another inert indicator such as 

30 helium. By choosing indicators that are not present m me body, me problem of accourmngfer 
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residual indicator is eliminated. That is, one may account for all of the indicator injected as it 

is released and exhaled. 

The above and still further objects, features and advantages of the present invention will 
become apparent upon consideration of the following definitions, descriptions and descriptive 
5 figures of specific embodiments thereof wherein like reference numerals in the various figures 
are utilized to designate like components. While these descriptions go into specific details of the 
invention, it should be understood that variations may and do exist and would be apparent to 
those skilled in the art based on the descriptions herein. 

RKTFF DESCTTPTTON OF T HP DP AWINGS 
10 Fig. lis a graph showing Ihe linearity and accuracy of me gas andvza o^ 

output monitor of the present invention over a range of O*, N*. and CO* and concentrations. 

Fig. 2 is a graph depicting typical measurements of rebreathed O*. N*, and CO, measured 
by the gas analyzer of the cardiac output monitor of the present invention. 
Fig. 3 is a graph showing the Hamilton respiratory flowmeter 
15 as a function of density for Oj, N 2 , and CO* 

Fig. 4 is perspective view of an experimental (breadboard) metabolic sensor used 
measure cardiac output in accordance with the present invention. 

Fig. 5 is a graph showing typical measured respiration parameters, including 
instantaneous respired flow rate and O2 and CO2 concentrations. 
20 F ig.6isagraphfflustrati^ 

and total volumes consumed and produced for a single bream. 

Fig. 7 is a graph showing the single breath CO* versm O2 concentrations fw me aWeol^ 

plateau region. 

Fig. 8 is a graph fflustrating a typical alveolar CO* concentration versus instantaneous 

25 respiratory exchange ratio. 

Fig. 9 is a graph showing cardiac output versus O2 consumption resulting from a variety 

of research studies. 

Fig. 10 is a graph showing expired CO* fraction versus expired volume mat is useful for 
determining deadspaos. 
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Fig. 1 1 is a flowchart illustrating an iterative technique fox determining cardiac output, 
shunts and deadspace in accordance with an exemplary embodiment of the present invention. 

DESCRIPTION OF THE PREFERRED E MBODIMENTS 
The following detailed explanations of Figures 1-1 1 and of the preferred embodiments 
5 reveal the method and apparatus of the present invention. 

The main components of the cardiac output monitor of the present invention are: alow- 
cost respiratory gas analyzer, a respiratory gas flowmeter, and the appropriate numerical 
algorithms necessary to make the calculations of cardiac output and the physiological corrections. 
The multiple medical gas respiratory gas analyzer, examples of which are disclosed in the 
10 aforementioned Drzewiecki patent applications, has the capability to quantify gas concentrations, 
including inhaled and end-tidal concentrations (approximating arterial and mixed venous blood 
partial pressures), of any constituent of respiratory gas mixtures of a known number of possible 
constituents, in real time on a breath-by-breath basis. A respiratory flowmeter, described 
hereinbelow, accurately determines the volumetric and mass flow rates of any gas/gasses as 
15 calculated from me product of measured total respiratory flow and the measured volumetric 
concentration. With this flowmeter, the resulting inhaled/exhaled volumes of the respiratory gas 
mixture are quantified in real time on a breath-by-breath basis. 

The gas analyzer of the cardiac output monitoring system of the present invention can 
determine in real time the individual concentrations of fluid constituents in a mixture ofN fluids 
20 by measuring independent properties of the mixture. In particular, N equations that, from first 
principles, relate the individual fluid concentrations to measured properties of the mixture, are 
solved for the N unknown individual concentrations of the fluids in the mixture. N-l properties 
ofthe mixture are measured by N-l sensors, which from cost considerations are preferably fliridic 
sensors, but may be any other technology devices, and N-l ofthe N equations are formed from 
25 me determined properties. The Nth equation is the constitutive equation which requires that the 
sum ofthe unknown concentrations of the N known constituents be equal to unity. 

For example, as described in greater detail in the aforementioned Drzewiecki patent 
applications, the individual concentrations of four gasses in a mixture of four known gasses can 
be determined by measuring the ambient pressure, temperature and flow rate ofthe sample flow 
30 ofthe mixture, the subsequent pressure drop ofthe mixture sample flow across a capillary and 
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an orifice whid.rn.yb* the supply nozzle of the flowmeter oscillaor. and find* ft. 
acoustic velocity^ nurture. Tnesa^eflowratecanbeme^bypass^^flow 
toough a fluidic feedback oscillator and measuring the output fteqnency period wtoeh » 

four unknown concenuations of the four known gasses ate dete^ by solving in real tm» 

fcurmdependerfe^uonsC^ane^ 

^uationretatrngmixtureviscosiryto^ 

heat to concentrations, and the constitutive equation). ■ 
10 Preferably, the oscUbtor flowmeter, sonic oscillator and the caputary are formed a. a 

mMs urements can be performed. Any one of the escalator nozzles can serve ..ft. ortfic* 
^dtah-»ht*»«-ftT.«— •«•»«. Tnedisposablesensormc^>seonnec*d 

^pBfiersusedtomeasureftechar^sticsofftenux^ 

Hne for drawing a sample. ... . , 

Advantageously, tow cost, fluidic so.sors measure flow, densny.v^ty airspeed 

ofMu odmgasmi*»re* 1^^**^^ 
20 transducers, lowest 

'e^ectriciummicrophonespr^ 

of the present invention requires no user calibration or maintenance and may be integral into 
^ monitoring systems For example, the gas analyzer can be added atongm. same flaw 
oam as other sensors or cm be added ma separate flow paft. 

« AKtougtaflwdiesensorsaremtf^ 

canbeimplementedv^ofterwesofsenson, For example, piezo^ecrtcaU^drrven suxft- 

acousucw^CSA^devicesh^beo^us^ 

devices can densfty. and electro-chemical devices can measure viscosity. Dependmg on theu- 
.»« and accuracy advantoges, these devices rnay be advantageously used m ptace «f 

30 fluidic sensors. , , . 

One of the important advantages of me g^ analyzer ofpreaent invents Btt»abm1y to 
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simultaneously determine the individual concentrations of N gasses in a mixture of N known 
gasses by using inexpensive sensors to measure properties of the mixture as a whole and by 
solving Nmdependentequauonsrel^ Tne number of gasses 

whose individual concentrations can bedetennined canbemCTeasedby mcorpo^ 
5 gas analyzer addiuonal sensors that measure additior^ 

a whole. If additional properties of me inixture can be independently measured by any means 
and related to unknown concentrations, concentrations of additional gasses can be determined 
In general, if N-l independent properties of the mixture of gasses can be measured, then N 
equations can be developed and solved for N gas concentrations (the Nth equation being the 

10 constitutive equation). 

Operation of the gas analyzer of the cardiac output monitor of the present invention has 
been experimentally verified. On a dry gas basis, a patient exhales a rmxture of O* COjimd 
enriched N 2 (a fixed rnix of N* CH* Ar, and trace gasses). A very low sidestream flow (~40 
ml/nun) sensor comprising a capillary viscometer and an orifice densitometer has demonstrated 

15 analysis accuracy of better than ^vol^ 
andQ2. Fig. 1 shows toe accuracy and 

twelve randomly selected known-assay mixtures of Cfe, O)* and N 2 . Analysis accuracy and 
resolution for gasses with significantly different physical properties, such as the volatile 
inhalation anesthetic agents Halothane and Isoflurane, in rnixtures of O* CO* is ±0.05vol% with 

20 a resolution of less than 200 ppm. Response is real time Cr90<200ms). Typical responsivny is 
evidenced by Ihe capability to measure breaflungin realtime. T^icalofthisismerecordofa 
rebreathing procedure shown in Fig. 2. Starting with pure O* me Oz decreases as CQ, 
concentration increases, and N 2 increases as the tissues excrete N 2 during denitrogenization. 
(Note that COj does not equflibrate probably becauserecircula^ 

25 ExperimenMresmtsinmcatematmegasa^ 

agents with resolution and accuracy comparable to those of IR devices, assays of O* mat are 
significantly more accurate than that of fael cell, paramagnetic and Oark electrode C>2 
and assays of N 2 that are superior to Raman systems. Also, with the exception of mass 
spectroscopy and Raman scattering, this gas analyzer is the first to offer real time quantification 

30 ofmtrogen, which can be* ^ 

circuit leaks and disconnections, as well as air emboli. 
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In addition to gas concentrations, which can be calculated using the aforementioned gas 
analyzer, a determination of cardiac output from respired gasses using a Fick-based technique 
typically involves a determination of the respired flow rate. For example, certain Fick techniques 
xequire the accurate quantification of both oxygen consumption and carbon dioxide production. 
5 The rate of consumption, or production, is the product of the instantaneous concentration and 
the flow rate. Measurement of respired flow rate (pneumotachometry) is conventionally made 
with one of a variety of devices: turbine meters, rotameters, fixed or variable orifices, capillaries, 
hot wire/film anemometers, ultrasonic/acoustic transit time sensors, etc With the exception of 
ultrasonic/acoustic devices, the accuracy of the flow through a fixed (nozzle, venturi) of variable 
10 (rotameter, flap) orifice depends on an intimate knowledge of, primarily, density as noted in the 
Bernoulli orifice equation, 

AP-pOWA 2 ) (1) 

where &P is the pressure drop, p is the density, Q is the volumetric flow, ca is the discharge 
coefficient (which typically is viscosity dependent) and A is the cross-sectional area. (A similar 

15 density-dependent equation can be derived for turbine meters where me density dependence 
comes from the conversion of the fluid kinetic energy to motion of the moving part/vane.) 

Since the gas analyzer of the present invention measures density and viscosity of the 
respired gas, a very low cost flowmeter, such as a bi-directional Hamilton variable area orifice 
device can be used to accurately compute Q because density is inherently known. Thepressure- 

20 flow relationship for this device for the three constituent gasses in air (nitrogen, oxygen, and 
carbon dioxide), is shown in Fig. 3 which clearly demonstrates the device's density dependence 
(e.g^ high density CO2 has lowest flow at the same pressure drop) and the error incurred if 
density is unknown. The pressure-flow relationship for this device, where the im» is a function 
of the pressure drop via the displacement of a wedge-shaped flap can be shown to be: 



25 



Q 2 - [L^pVp] [AP 3 - AP'MAf W] (?) 

where L is the characteristic dimension of the flap, Af is the area of the moving flap andk is the 
effective spring constant of the canmevered flap, and p , again, is the density. 

Given the accurate determination of volumetric flow, the product of the flow, Q, with the 
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individual gas concentrations during the exhaled breath gives CO* production and Oi 
consumption, the ratio of which is respiratory quotient, R, a term that is critical in me 
determination of cardiac output, as detailed hereinbelow. 

Fig. 4 illustrates an experimental cardiac output monitor 10 constructed to demonstrate 
5 the feasibility of readily acquiring the data required to implement the methodologies of me 
present invention. Cardiac output monitor 10 includes a face mask 12 for receiving respired 
gasses from a subject A respiratory flowmeter 14, coupled to the face mask via a humivent 16, 
measures the volumetric respiratory flow. A sample pott 18 at the output of flowmeter 16 
provides a sample flowto the gas analyzer. The gas analyzer includes a gas sensor 20 contaimng, 

10 for example, meaforemenuon^^ 

a transducer module 22. A vacuum line 24 draws the sample gas through the gas sensor 20 of 

the gas analyzer. The transducer module 22 and flowmeter 14 provide measurements to a 

processor (not shown) which determines cardiac outout and pulmonary 

with Fick techniques implemented in software. 
15 By applying a mathematical algorithm that combines and integrates concentration and 

volume over time, both the rate and amount of uptake (gain) or release (loss) of any constituent 

(indicator) can be determined in quasi-real time, that is, once every several breath cycles as 

opposed to real time which implies instantaneous^ contmuous reading 

blood flow pumped by the heart) can thus be calculated using the well-established and accepted 
20 Fick Principle in real time using only breath-by-breath concentration information. It should be 
understood that a number of different non-invasive techniques for measuring cardiac output from 
respiratory gasses based on the Fick Principle have been proposed over the years. The system 
of the present invention can employ any one or combination of these Fick-based techniques, 
including modifications and improvements mereto, to deternmie cardiac outn^ 

Because of the ability to measure the rate and amount of uptake as well as the end-tidal 
partial pressures of unique indicators such as nitrogen, nitrous oxide, anesthetic agents, and other 
inert gasses in real time, and/or concentrations of respired gasses such as oxygen and carbon 
dioxide, it has been possible to develop mis unique noninvasive cardiac output monitor using the 
Fick Principle. While conventional gas analysis technology, such as IR spectroscopy, could be 
30 used to measure nitrous oxide and anesthetic agent concentrations, it is the unique ability to 
inexpensively measure nitrogen and other inert gasses such as helium that offers both a non- 



25 
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invasive but also a non-afFecting, non-toxic approach. Furthermore, it is the ability to measure 
the density and viscosity properties of the inspired and respired gas mixture (e.g., with other 
technology one would have to measure the concentrations of all, not just some, constituents 
sunultaneously) that allows for the measurement of the inspired and respired volumetric flbwrate 

5 independent of the properties. By combining these gas analysis measurement techniques with 
accuratemodels and simulations of uptake and distribution of the gasses in the body, the current 
concept can be further refined to consider pulmonary shunts and deadspace as well as other 
individual-specific physiological and disease states. 

In the invasive appHcation of the Fick Principle technique classically referred to as me 

10 indicator dilution technique, a known amount of an easily detectable, substance (indicator) is 
injected into the bloodstream and its presence is monitored a short, known, distance downstream. 
When the indicator and blood are completely mixed, the concentration-time relationship of me 
indicator provides information to deterrmne the blood flow rate, which, by definition, is the 

cardiac output. 
15 jn accordance wim one emb<>dimeiuo 

the Fick technique is used and is suited to analysis of respiratory gasses with the gas analyzer of 

the present invention. Specifically, rather than directly (invasively) injecting an indicator into 

the bloodstream, a known volume/mass bolus of an indicator gas (e.g. nitrogen, nitrous oxide, 

sevoflurane, desflurane, or helium) is introduced mto me in^^ 
20 patientbreatbingcircuit Tbe blood will take up this gas via respiratory transfer in the alveoK of 

fcehmg* In order to rninirnize the time it takes for the indicator to fill the alveolar space and 

to reduce the time required for the diffusion of me mdi^ 

ventilated at a redetermined, suffident, breathing rate. Tnebolu^ 

during one or more quick inspired breaths. Hie blood will take up only a portion of the bolus of 

25 respiratory mdi^ 

directly into the bloodstream. The transfer of gas to the blood occurs because there is a 
difference in the partial pressure of the gas in the alveoli and in the blood stream. The alveolar 
partial pressure is higher than that in the blood, there being no (or a different amount of) gas m 
the blood. In the first inspired breath/breaths, the indicator enters the lungs. Subsequent 

30 inspired breath wffl be rreeofmtotorbu^ 

start to dilute and evacuate the lungs. Within the first four or five exhaled breaths, the indicator 
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gas that has not been taken up in the blood fiom the lungs will be completely (to over ninety-nine 
percent) removed. 

The amount of indicator exhaled in these initial breaths is measured using the property 
detecting, multiple gas analyzer in conjunction with a respiratory flowmeter inline with the 
breathing circuit The amount of indicator gas taken up by the blood is determined by subtracting 
me measured exhaled amount from the known input bolus. It is valid to perform this 
computation and measurement in the four or five breaths it takes to purge the lungs of the 
indicator gas, because there is no release of the indicator from the blood since the blood has not 
yet returned to the lungs after circulating through the body. The indicator gas is chosen to be one 
that has a low solubility in the tissues and so remains primarily in the blood. Nitrogen is such 
a gas. When the blood carrying the indicator gas returns to the lungs, there now exists a reverse 
difference in partial pressure between the blood and the alveoli. This time, the partial pressure 
in the blood is higher because there is no indicator in the alveoli. The end-tidal values (which 
approximate the arterial blood gas partial pressure) of indicator gas in each exhaled breath are 
measured and the amount and time history of the gas released from the bloodstream after it has 
circulated through me cardiovascular system is simu^^ The total volume of exhaled 

indicator and the time it took to be released from the bloodstream, the area under the exhaled 
indicator gas-time curve, compared with the amount of indicator taken up, can be used to 
calculate the cardiac output using the Stewart-Hamilton relationship if no shunt or deadspace 
effects are present The Stewart-Hamilton relationship is discussed and derived hereinbelow. 

Indicator dilution techniques can be subdivided according to the method of application. 
Two commonly used techniques are continuous infusion and bolus injection. The continuous 
mfusion method was first described by Stewart A major disadvantage of mis technique is that, 
in the closed circulatory system, the indicator eventually saturates the blood stream and any time 
history information is lost The injection of a bolus of indicator, on the other hand, does not 
result in such saturation, and, although lower concentrations must be resolved, is used more often 
in clinical practice and is the technique that is a basis for a preferred embodiment of the present 
invention. 

The principle of both methods can be illustrated by means of a simplified flow system in 
which a constant flow of fluid with a flow rate Q 0 and uniform velocity profile is assumed for a 
system with a single inlet, injection site, and outlet, detection site. Ideal, 0*., complete) cross- 
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sectional mixing is assumed throughout the volume of the system between injection and detection 
site. 

For continuous infusion of indicator gas with injection rate oj into a steady blood flowrate 
Qo> ^ volumetric concentration measured, Co, after complete mixing has occurred, is constant 
and determined by the ratio of volumetric flowrate of indicator to flowrate of blood, 

Co=qyQ e < 3 > 
or, after rearranging and solving for the flowrate (e.&, cardiac output), 

Qo-oyco. (4) 

When a bolus of indicator is injected into the blood streain, me relation between flbw^ 
and concentration is obtained by using me condition of conservation of indicator. In an 
infinitesimally small time interval dt, the volume of indicator dV lP assing out of me ouflet equals 
the concentration of indicator at that time, c<t), multiplied by the volume dV b of the blood 



dVi = Ci(t)dVb 



<5) 



15 The blood volume dV„ equals Ihe blood flow rate, Qo(t), multipUed by me time interval dt, 
leading to 

dVi-Qo(t)c(t)dt <«> 

With use of the condition of conservation of mass, the total volume of indicator injected will pass 
me detection site, yielding: 



20 



V, = jQo(t)ci(t)dt 
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For constant flow, the flow rate can be separated from the integral and solved for so that, 



<«) 



Jci(t)dt 



0 



This equation is referred to as the Stewart-Hamilton equation. From this equation, the 
flow rate can be determined by the quotient of the total amount of indicator gas injected and the 
area under the measured indicator dilution curve. The ctf) curve has an asymmetrical shape, 
which is caused by the continuation of the mixing process between indicator and transport 
medium, even when the first portion of indicator has already passed the detection site. 

In the non-invasive dilution-indicator technique of the present invention, the indicator is 
a known gas constituent of a respiratory gas mixture. Preferred indicators are: nitrogen, nitrous 
oxide, helium and a relatively insoluble volatile anesthetic (sevoflurane, desflurane). At some 
initial time zero, a bolus of indicator gas is injected into the respiratory breathing gas and 
ventUatedmto toepatierf^ Following injection, 

the ventilatory cycle is continued and the exhaled gas is anaryzed and amount of me m^ 
exhaled during each breath is quantified. From the amount of gas measured in me exhaled 
volume, the content of indicator in the blood are calculated using Henry's Law. The measured 
concentration of indicator as a function of time elapsed is integrated over time to determine the 
total amount of indicator exhaled during the time interval. By incorporating mis data (amount 
of indicator taken up by the blood and the area under the exhaled indicator time curve) into 
equation (8), the cardiac output is calculated. This technique is analogous to both dye and 
mermaldnutionmemods,butbasmeadvanteg^ ■ 

The assumptions used in the derivation of the Stewart-Hamilton equation for indicators 
include: (1) an open flow system (i.e., the indicator passes through only once); (2) no loss of 
indicator (i.e., no uptake in tissues); (3) complete cross-sectional mixing resulting in a uniform 
concentration distribution; and (4) constant flow (no transients). In clinical situations these 
conditions are rarely fulfilled when using direct mjection of an mmcatorsubsto 
stream. Since 1he circulatory system is necessaruy closed, recfc^ 

Depending on the solubility of the indicator, loss of indicator cannot always be avoided. 
Depending on me injection and detection sites, adequate mixing may not always be achieved. 
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And, most importantly, constant flow does not exist due to the pulsating cardiac action and 
effects of spontaneous or mechanical ventilation. 

These factors are ove^mem the present invennon. By introducing into the blood stream 
a relatively small volume of indicator gas that is clear of indicator over a short period (i*, less 
tban the recirculation time), recirculation does not come into play. By using an insoluble gas,- 
either nitrogen or very small amounts of relatively insoluble agents such as helium, nitrous oxide 
or sevoflurane or desflurane, the loss of indicator is extremely small or negligible and again any 
effects of recirculation are not significant Since nitrogen is relatively insoluble compared to the 
other indicators, and the tissue compartments and venous blood are saturated, at equilibrium and 
constant, changes in arterial blood nitrogen content reflect directly the changes in the arterio- 
venous difference. By introducing the indicator into the respiratory gas during inhalation, rnixing 
is not an issue, because the transfer through the microscopic capillaries in toe alveoli ensures 
umformityin the blood stream. Finally, by averaging at least three esthnates over the ventilatory 
cycle, the effects of variable flow are minimized. 

Although it is possible to accurately measure the variables described above, additional 
consideration is required to account for the presence of pulmonary shunts and deadspace in the 
lungs. Shunts are regions where blood vessels (capillaries) bypass the alveoli and as a result do 
not permit blood to come into contact with the indicator gas. Thus, in the presence of shunts, any 
cardiac output that is computed by the above-described procedure, while accurately computing 
the blood flow to the alveoli, could mderestimate the actual blood flowrate because the shunt 
flow is not determined. Deadspace, on the other hand, is a region of the lungs where no transfer 
of indicator (or oxygen, etc) takes place. This relates to what is known as lung capacity. Alveoli 
that do not contribute to blood gas transfer serve only to dilute the end-tidal value of me 
indicator. That is, when indicator is given up by the blood in working alveoli, gas with the 
indicator mixes with gas without the indicator from the deadspace. Therefore, in the presence 
of deadspace, any computation of cardiac cutout by the above-described procedure would 
overestimate blood flow because of lower end-tidal readings which suggest shorter residence 
time in the lungs of the blood and, hence, higher blood flow or cardiac cutout The inability to 
compensate adequately for these two factors (i.e., shunts and deadspace) is a major shortcoming 
of indirect Fick techniques and probably the reason why such techniques have not been 
successful in the past 
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In accordance with one embodiment of the present invention, in order to overcome mis 
serious deficiency and at the same time provide a method for quantifying shunts and deadspace, 
the measured time-history of exhaled indicator is compared to results for the same time-history 
from a detailed physiological model of Ihe human cardio-vascular system that contains all me 
5 necessary parameters to describe the uptake, dista^ution and release of respired gasses. The Gas 
Uptake Distribution (GUS) model, developed by one of the present inventors, is a computer 
model that meets this requirement The GUS model simulates respiratory and anesthetic gas 
uptake and distribution using a simplified, but functionally accurate, anatomical and 
phyaologicalmodelofthehumanbody. The simulation utilizes the patient parameters of weight 

10 and percent body fit, in conjunction with a parametric description of patient physiology <e^, Cfe 
consumption, pulmonary shunt, systemic shunt, cardiac output, deadspace), ventilation mode 
(controlled, spontaneous, minute ventilation) and gas delivery (air, Oj, N 2 6, anesthetics). Real 
time predictions of the concentrations, and uptakes or losses of any one of nine respiratory and 
anesthetic gasses in twelve compartments can be observed. The twelve compartments consist 

15 ofeleven tissue compartments and one patient breaming circuit of ^ 

or mechanical ventilator. While the exemplary GUS model described herein involves twelve 
compartments, it will be understood that the GUS model can involve any number of suitable 
compartments sufficient to accurately model respiratory and anesthetic gas uptake and 
distribution in order to quantify the effects of deadspace and shunts. 

20 The GUS model is integrated with the gas analyzer system microprocessor and is used 

to predict, in real time, the rate of exchange of an indicator between the alveolar gas and blood, 
the end tidal partial pressure of the indicator during normal and non-uniform ventilation. The 
simulation uses inputs based upon known patient data to make an imtid prediction of alveolar 
partial pressure of me indicator. This initial approximation is then modified to predict the 

25 measured rates of uptake (gain) and release (loss) and end-tidal partial pressures of the indicator. 
Because the respiratory gas analyzer, in conjunction with a respiratory flowmeter (spirometer), 
can actually measure metabolic rates and O2 consumption, and the ventilation parameters are 
known inputs, the only unknowns are the physiological parameters of shunts, deadspace and 
cardiac output These are manipulated in real time using fuzzy logic as me si^ 

30 to develop an output matching the measured parameters. The resulting values of shunts, 
deadspace and cardiac output represent an accurate estimate of the actual values. Once these 
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values have converged at some value of cardiac output that incorporates shunts and deadspace, 
that value of cardiac output is reported as the output of the monitor. 

The present invention offers significant advantages over olher technologies that utilize 
the Fick Principle for detennining cardiac output As previously explained, the system of the 
present invention has the capability to determine cardiac output using any of the many recognized 
Fick techniques By implementing the present invention, the clinician is not limited to a single 
application of the Fick Principle but is offered more flexibility by selection of appropriate 
techniques, methodologies, indicators, and algorithms for me calculation of cardiac output that 
fit the clinical situation. For example, currently available methods are limited to only one 
indicator (carbon dioxide or oxygen) at a time (e.g., the NICO system). With the present 
invention, any respiratory or anesthetic gas may be used as an indicator and more than one gas 
can be used simultaneously. This enables the clinician to select the appropriate indicator for the 
clinical situation. 

The determination of the rate of uptake (gain) or release (loss) of a substance (indicator) 
in the alveolar gas is determined by multiplying the concentration of the substance (indicator) at 
the subject's airway by the airway flow and integrating the product over the respiratory cycle. 
Using this procedure the inspired and expired volume (mass) of the substance (indicator) is 
deternimedon abreath-by-breathbasis. The gain/loss of the substance (indicator) can then be 
calculated as the difference between the inspired and expired substance volumes divided by the 
duration of each breath. This can be expressed as: 

V^^a/rd/V^cXt^-ll/fTB-WVKt)^^ {integrals) (9) 
where VKt) = airway flow volumeper unit tune 
Cj(t) = airway substance concentration 
Ti — end of inspiration (time) 
Tb *» end of expiration (time) 

The fiaction of the substance (indicator) at the patient's airway, c,<t), is multiplied by the airway 
flow, Vi(t), and the product is integrated over the respiratory cycle. Using this procedure, me 
inspired and expired volume of the constituent (indicator) is determined on a breath-by-bream 
basis. The uptake (gain) or production (loss) can then be calculated as the difference between 
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the inspired and expired volumes of the substance (indicator) divided by Ihe duration of each 
breath. These values are then inserted directly into the Fick Equation for calculating cardiac 
output 

The Fick Principle described hereinabove was first articulated in 1 870 by Adolf Fick and 
5 is easily derived from the basic laws of transport phenomena with the appUcation of the 
Conservation of Mass to blood flowing through an organ system. The Tesult of this derivation 
simply states that the flow of blood through an organ is equal to the total uptake or release of any 
substance (indicator) by the organ divided by the arteriovenous concentration difference of the 
substance (indicator). Another way of expressing the Fick principle is that the size of a stream 

10 may be readily calculated if we know the amount of substance (indicator) that enters or leaves 
the stream and the concentration difference resulting from such entry or removal. Both of these 
statements are true whether the substance (indicator) is a respiratory gas such as oxygen, carbon 
dioxide, nitrous oxide, or an intravenously injected substance such as cold normal saline. 

The application of the Fick Principle to the measurement of pulmonary capillary blood 

15 flow enables the clinician to approximate cardiac outoutdi^^ 

levels of pulmonary shunting ofblood. A mathematical application of the Conservation of Mass 
to a cardiopulmonary gas exchange model yields an expression for direct and indirect Fick 
techniques that has become known as the Fide Equation: 

Q»=N7(Ca-Cv) (10) 
20 where: Q»= pulmonary blood flow (cardiac outom) 0 Mood^mm) 

jj» «_ rate of gain of substance (indicator) from alveolar gas (ml/min) 
Ca - concentration of substance (indicator) in arterM blood (ml/1 blood) 
Cv = concentration of substance (indicator) in mixed venous blood (ml/1 blood) 
In order to apply the Fick Equation to determine pumionary blood flow, techniques mat 
25 can either directly measure or indirectly estimate the rate of exchange of a substance (indicator) 
between the blood and alveolus must be readily available. The same is true for deterrnining the 
concentrations of the substance (indicator) in both the mixed venous and arterial blood. By 
incorporating this data into the Fick Equation, it is possible to calculate the pulmonary blood 
flow and thus cardiac output. 
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Numerous methods have been developed to determine each of the variables required for 
the calculation. These methods vary in computational complexity, equipment utilized, and the 
supervision required. Some methods involve Evasive techniques (direct Fick) while others 
utilize various empirical correlations between noninvasive measurements and the variables of 
5 interest (indirect Fick). Various indicators including oxygen, carbon dioxide, nitrous oxide, and 
anesthetic gases have been used in these respiratory Fick techniques. 

The use of direct (invasive) Fick techniques provides the most accurate prior known 
measurement of cardiac output Unfortunately, direct techniques require sampling of both mixed 
venous and arterial blood in order to determine the amount of indicator in each. This requires 
10 the insertion of needles and catheters into the systemic circulation and heart in order to obtain 
the specimens required for analysis. 

Indirect (noninvasive) Fick techniques rely on measurements of alveolar gas 
concentrations to estimate the amount of each indicator in arterial and mixed venous blood. 
Indirect (noninvasive) determinations of arterial substance (mdicator) cont^^ 
15 based on the partial pressure of the gas sampled at the patient's mouth at me end of a tidal 
expiration* 

For the indirect Fick techniques, two assumptions are made. The first assumption is that 
the partial pressure of the gas sampled at the end of a tidal expiration (PET) provides a good 
representation of the alveolar gas. Secondly, the alveolar gas is in equilibrium with the substance 
20 (indicator) in the arterial blood leaving the lungs. If these assumptions are valid, PiET can be 
used to estimate the arterial partial pressure of the substance (indicator). However, it is important 
to emphasize that by mcorporatm^ 

output is limited during conditions a large pulmowuy shunt and/or a la^ 
fraction. When these conditions exist, appropriate corrections to the estimation of cardiac output 
25 must be made. 

To correct for the potential effects of a large pulmonary shunt and/or alveolar deadspace, 
a somewhat more detailed model of the cardiopulmonary system than me one used for the 
derivation of the Fick Equation is necessary. This model must allow for pulmonary shunting of 
blood flow and include an alveolar deadspace compartment A pulmonary shunt represents 
30 pulmonary blood flow that perfuses regions of me lungs that are not being vermlated. An 
alveolar deadspace compartment is defined as a region of the lung mat is being ventilated but is 



PCT/USOO/01465 

26 

Application of the Fick principle for the indicator at the blood-lung interface 

QT-Q , s = QT(l.{Q , s/QT}) = N/(Cv-CaI) (11) 
where QT = Total pulmonary blood flow (ml blood/min) 

Q f s = Blood flow that does not xmdergo gas exchange - shunt (ml blood/min) 
N = Loss of indicator to alveolar gas (ml/min) 

Cv = Concentration of indicator in mixed venous blood (ml/ml blood) 
Cal = Concentration of indicator that has undergone gas exchange of indicator 
with alveolar gas (ml/ml blood) 

Similarly, a mass balance of the indicator at the subject's airway, assuming that the 
alveolar deadspace empties in parallel, results in: 

Cm*VT==CA(VT-VD) + CD*VD (12) 
where CA = concentration of indicator in alveolar gas (ml/ml gas) 

Cm ■« Concentration of indicator measured at the subject's airway (ml/ml gas) . 
CD = Concentration of indicator measured in deadspace compartment (ml/ml gas) 
VT = Tidal volume (ml) 
VD « Alveolar deadspace (ml) . 

Solving for CA yields: 

CA = (Cm*VT- CD*VD)/(VT-VD) (13) 

With the utilization of the various measurement techniques to determine the appropriate 
variables for these equations, cardiac output can be calculated with correction for any effects of 
shunt and deadspace incorporated into the final value. 

The measurement of the variables used to calculate cardiac output is accomplished by 
utilizing three common Fick techniques: direct Fick, indirect Fick, and indicator-dilution. The 
variables to be measured are the uptake/release (gain/loss) and the arteriovenous concentration 
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difference of the indicator. The present invention enables the clinician to measure these variables 
in numerous clinical situations. This approach offers significant advantages over other methods 
utilizing the Fick Principle for determining cardiac output 

As discussed above, in order to apply the Fick Equation to determine pulmonary blood 
5 flow (cardiac output), the rate of exchange of a substance Ondicator) between the alveolar gas- 
and blood and the arteriovenous difference of the substance (indicator) must be measured. Each 
of these variables is detennined from the measurement of two other variables. The rate of 
exchange of the indicator can be determined by measuring the amount of indicator (constituent) 
that is present in both the inhaled and exhaled breath. The difference in amount between inspired 
10 and expired constituent is the net gain or loss of the indicator that is exchange between the 
alveolus and the blood. The arteriovenous difference of the indicator is determined by 
measuring/estimating the amount of indicator (constituent) that is present in both the arterial and 
mixed venous blood. The difference in these two values is the arteriovenous difference in die 
blood. 

15 In accordance with the present invention, the determination of the rate of uptake (gain) 

or release (loss) of a substance (indicator) in the alveolar gas is determined by multiplying the 
concentration of the substance (indicator) at the subject's airway by the airway flow and 
integrating the product over the respiratory cycle. Using this procedure the inspired and expired 
volume (mass) of the substance (indicator) is determmed on a bream^y-bream basis. The 

20 gain/loss of the substance (indicator) can then be calculated as the difference between the 
inspired and expired substance volumes divided by the duration of each breath. This can be 
expressed as: 

V'(t)^ 1/TI// V'(t)F(t)dt- {1/(TE-TI)}// V'(t)F(t)dt {integrals) (14) 
where V*(t) = airway flow 
25 F(t) = airway substance concentration 

TI = end of inspiration (time) 
TE = end of expiration (time) 
A microcomputer is used to multiply the fraction of the substance (indicator) at the patient's 
airway, F(t), by the airway flow, V(t), and to integrate the product over the respiratory cycle. 
30 Using mis procedure, the inspired and expired volume of the constituent Ondicator) is determined 
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on a breath-by-breath basis. The uptake (gain) or production (loss) can then be calculated as fee 
difference between the inspired and expired volumes of the substance (indicator) divided by the 
duration of each breath. This value is then inserted directly into the Fick Equation for calculating 
cardiac output. 

5 Because of the increased capability and versatility of the present invention, at least three 

unique techniques for measuring cardiac output noninvasively using the Fick Principle are 
possible. The first is an adaptation of the Soluble Inert Gas Technique. The second is an 
adaptation of the indicator-dilution method. The third is an integrated predictive computer model 
technique incorporating the computer simulation (e.g., the aforementioned Gas Uptake 

10 Simulation (GUS) model). As previously described, the GUS model can also be used to make 
correction for the presence of increased shunt and deadspace. 

Several noninvasive respiratory applications of the Fick Principle are based upon the 
uptake and distribution of physiologically inert and soluble gases. When a soluble gas is inhaled, 
the partial pressure of the gas in the pulmonary capillary blood equilibrates with the partial 

15 pressure of the gas in the alveolar space. The volume of gas absorbed by the blood can be 
calculated directly from the partial pressure of the gas in the alveoli, the solubility of the gas in 
blood, and the quantity of blood that has equilibrated with gas. This can be expressed 
mathematically using Henry's Law as: 

Vgas = Sgas*PAgas*Q ( 15 > 

20 where: 

Vgas - volume of gas absorbed by the pulmonary capillary blood (ml) 
Sgas «= solubility coefficient of the gas in blood (ml gas/ml blood/mm Hg) 
PAgas = alveolar partial pressure of gas (mm Hg) 
Q = volume of blood equilibrated with alveolar gas (ml blood) 

25 Differentiating the expression with respect to time and rearranging terms result in 



q» t= V'gas /(Sgas*PAgas) 



(16) 
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Equation (1 6) represents a special case of the Fide Equation where the mixed venous gas 
concentration is zero. The absence of the mixed venous term is due to the fact that under normal 
conditions me indicators used are not present in blood. 

The use of a soluble inert gas as the indicator in the Fick Equation therefore eliminates 
the need to determine the mixed venous gas concentration provided that the mixed venous gas 
concentration is zero at the start of the measurement procedure. However, this zero concentration 
assumption requires that the measurement procedure be completed before the blood recirculation 

tunc* 

Primarily two gases have been used as indicators in inert gas techniques, nitrous oxide 
and acetylene. For the present invention, nitrous oxide and the newer volatile anesthetic agents 
(sevoflurane, desflurane) are preferred over acetylene. These gases have a relatively low 
solubility in lung tissue and higher solubility in blood. Since high concentrations of mese gasses 
in fhe blood are undesirable, low levels of inspired gas concentration must be used. This low 
level of inspired indicator concentration decreases the magnitude of both V'gas and PAgas in 
Equation (lo), ™ d therefore increases the sensitivity of the blood flow estimate to error in the 

determination of these parameters. 

A number of techniques have been developed to increase the accuracy in measuring bom 
Vgas and PAgas. Ttese techniques have been discussed an^^^ 
rebreathing techniques and single-breath techniques. 

Oneof me limitations of mese techniques is the assumption that the venous concentration 
ofmemdicatorgasiszero. This assumption not only limits the duration of the measurement of 
V'gas and PAgas due to venous recirculation, but also requires a waiting period between 
determinations to allow m^^ However, it is possible 

to overcome this limitation by using nitrogen as the indicator. Since nitrogen is relatively 
insoluble compared to the other indicators, and the tissue compartments and venous blood are 
saturated, at equilibrium and constant, changes in arterial blood nitrogen content would reflect 
directly the changes in the arteriovenous difference. 

The above-described embodiment is a major improvement over existing non-invasive 
techniques in that it provides means for compensating for the effects of shunts and deadspace. 
However, the fuzzy logic used in estimating the effects of shunts and deadspace may not resolve 
certain ambiguities that may arise. Another, more general embodiment of the present invention 
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involves separately estimating the deadspace (both physiolbgical and anatomical) so that only 
shunts need to be iterated in the physiological model software. In addition, in some cases, it may 
be inconvenient to provide indicator gasses that need to be tracked; accordingly, the second 
embodiment utilizes the unique instantaneous relationships between oxygen and carbon dioxide 
to provide much the same information. 

In "Deadspace and the Single Breath Test for Carbon Dioxide During Anesthesia and 
Artificial Ventilation,- Br. J. Anaesth (1984), 56, 109, incorporated herein by reference in its 
entirety, Fletcher et al. teach that the anatomical deadspace (i.e., the deadspace of the passages 
in the lungs that bring the ventilatory gasses to the alveoli) can be accurately measured by 
measuring the total volume of expired gasses before any change occurs from that inspired in 
measured concentration of either oxygen and carbon dioxide. Further, Fletcher teaches that me 
physiological deadspace (U., the measure of the alveoli that do not participate in transfer of 
gasses to the blood) can be estimated quite accurately by measuring the volume of expired gasses 
during alveolar release. Using relatively simple algorithms, mis deadspace is related to the slope 
15 of the concentration-volume trace as end-tidal values are approached. 

Kim et al. teach in the aforementioned paper that the value of the mixed venous carbon 
dioxide partial pressure is related to the value of the exhaled carbon dioxide concentration at the 
point when the ratio of the change of carbon dioxide concentration to the change in oxygen 
concentration (Le., the instantaneous respiratory quotient) is approximately 03. This ratio value 
is a function of physiological conditions as well as disease state; however, by applying me 
method of Fletcher et aL to account for deadspace, the variation maybe limited to shunts alone. 
Since cardiac output (not including the effects of shunts) is thus defined by the Fick equation to 
be me ratio of the carbon dioxide produced to 

dioxide partial pressures, this portion of the calculated cardiac output may be introduced into the 
25 physiological model software (e.g., GUS) and a shunt value iterated that results in the respired 
breath waveform as measured. In cases where it is reasonable to assume that shunts do not play 
a major part (e.g. when they are not present or remain constant during the period of momtoring), 
the data from the gas analyzer and respiratory flow meter provide cardiac output directly, with 
accuracy that is superior to existing rnodified-Fick techniques. A detailed example of a modified 
30 Kim single-breath technique is now provided. 

Given that the system of the present invention can inexpensively and reliably measure 
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real time, simultaneous Oa and CCh concentrations (and N 2 if needed) as well as consumption 
and production right at the patient breathing circuit, single-breath techniques can reliably be use 
to non-invasively measure cardiac output and provide continuous detennination of cardiac output 
on a breath-by-breath basis. 
5 As described in the aforementioned paper, Kim et al. developed a method of estimating 

•true" mixed-venous and arterial CCh partial pressures by the analysis of simultaneous exhaled 
O2 and CO2 partial pressures (Feoj and Fecoj) on the alveolar plateau of a single extended 
exhaled breath. Given either O2 consumption or COj production cardiac output can then be 
estimated using the Fick equation. Several studies have been undertaken since then that have 

10 attempted to validate and improve the technique, mostly from a data processing perspective. 
However, this methodology heretofore has not received wide acceptance primarily for three 
reasons. The first reason is the difficulty associated with obtaining simultaneous realtime Oj and 
CO2 data, putting the technique into a very high cost category. This limitation has been 
overcome in Ihe present invention with the use of the aforementioned real time multi-gas 

15 analyzer. The second reason is the extensive calculations that must be perfisrmed. This 
limitation has been overcome in the present invention by using state-of-the-art microprocessor 
technology. The third reason is the concern about the validity of the assumptions that were made 
in the original methodology. These concerns focus upon the potential effect of any body gas 
stores, alveolar deadspace and pulmonary shunts vis & vis disease state, me shape of the 

20 oxyhemoglobin and carbon dioxide dissociation curves, me necesdty for a relative steady ^s^ 
(also required theoretically for any Fick and single-injection indicator-dilution methods), and the 
existence of any ventilation/perfusion (V7Q) abnormalities (the effect of sequential emptying of 
the alveoli). Most of the assumptions made by the original investigators have subsequently been 
validated, but concerns have been expressed overusing the technique in patients with abnormal 

25 ventilation/perfusion ratios (VVQ). These concerns have generally been overcome by die 
significant number of studies that have categorically confirmed that there is sufficient 
information in the respired CCh data alone that correlates with cardiac output The former have 
been studies based on the initial efforts of Fletcher who identified a very accurate techmtme to 
account for deadspace, and have culminated with the very recent work of Arnold et aL who 

30 conducted detailed correlations of exhaled COj waveforms at various work/exercise levels to 
cardiac output These latter studies dearly confirm that me 
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elements relating to cardiac output 

In view of these studies, the present inventors have detennined that if the true 
instantaneous alveolar values of Cfe and CO2 that account for lung volume and physiological 
deadspace can be imputed/measured, the true instantaneous respirator exchange ratio can be 
5 determined directly for the alveoli, thereby allowing computation of the actual blood flow being 
passed through the lungs. Pulmonary and systemic shunting is accounted for with knowledge of 
the arterial oxygen saturation, which is routinely measured in a peripheral artery with a pulse 
oximeter. According to the methodology of the present invention, the single-breath technique 
is augmented to account for V/Q abnormalities. This is accomplish by modeling the physiology 
10 and, therefore, having a virtual patient on the side to compare with. The inputs of deadspace and 
cardiac output to the physiological model are iterated untfl me t^odel converges on the measured 
metabolic outputs. 

Using the cardiac output monitor 10 shown in Fig. 4, experimental data was collected and 
the computations required to determine cardiac output were performed, as shown in Fig. 5-8. 
15 Fig. 5 shows the typically measured respiration parameters, ir^ired and exhaled Cb and^^ 
concentrations and instantaneous resph^ flowrate, for several breathe Fro^ 
parameters are computed. Fig. 6 illustrates the rates of O2 and CO* consumption and production, 
their ratio (respiratory quotient), and total volumes consumed and produ 

Fig. 7 shows the single breath Fcoa plotted against F02 for the alveolar plateau region. This data 
20 is processed to obtain a relationship between alveolar CO2 concentration and instantaneous 
exchar^erau^^ Note 1hat me data spans the range 

of exchange ratio between approximately 02 and 0.8 so that only a minor extrapolation to the 
steady state respiratory exchange ratio (Rss) is required. Rss is calculated by integrating the 
product of the concentration and the flowrate. For this case, Rss - 0.889; however, it shouldbe 
25 noted that it also equals the ratio of the CQ* production to O2 consumption rates shown in Fig. 
6. Using the mixed venous and arterial values estimated from the points when R - 032 and R 
= 0.889, an arteriovenous difference of 0.9% (6.84 mmHg) is obtained. The minute oxygen 
consumption fiom Fig. 6 is 75.9ml m^^ This results ma cardiac 

outoutforthatparticularbream of approximately I0.n 
30 active state of the subject, and falls within expected ranges as. shown in the data of numerous 
other researchers in Fig. 9. 
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The diagnostic capability of this methodology is considerably greater than merely 
providing cardiac output Metabolic rate has been computed and can be reported to the 
user/clinician. Deadspace and pulmonary shunts can be quantified; information that would be 
of significant interest to the respiratory therapist as well as thepulmonologisL 
5 In Kim's technique, the lungs act as an aerotonometer that measures the partial pressure 

of the mixed venous blood, and hence gas content The essentid feature of Kim's method is me 
application ofthe Haldane principle to determine the mixed venous and arterial CO2 partial 
pressures (PVC02 and P A co2) from simultaneous measurements of COj and Q* concentrations from 
a single exhalation that is long enough to ensure that alveolar gasses have equilibrated and that 

10 enough data points can be collected on the alveolar plateau (where the exhaled gasses are 
representative of the alveolar gasses) for a statistically meaningful regression analysis to be 
performed. Henry's Law of Solubility and the dissociation curves relate areolar partial pressures 
and the content in the blood for CQ* and Q2. Using the cardiac outputmonitor of the present 
invention, with sample rates of 30 - 200 per second, data taken under normal respm 

15 8-12 breaths per minute yield satisfactory results; thus, patient cooperation is not required. 
Assuming that mixed-venous (Xh tension Pvco* remains constant during this time (blood 
normally does not recirculate in the period of one breath), men arterial CQj tension (Pace*) will 
rise and eventually, if the breath is long enough, will approach Pvcoa. Haldane observed that me 
arterial value is equal to the mixed venous value When the instantaneous respiratory exchange 

20 ratio, R - V002/V02 = 032. A detailed discussion of data supporting this observation is 
presented by Kim et aL and will not be repeated here. 

Kim, however, made the assumption 1hat the exhaled concentrations were equal to the 
alveolar values. Clearly, this would hold true only if there were no alveolar deadspace. The 
relationship of PCO2 with R should, of course, be correct for the air exhaled from a fuDy 

25 ventilated and perfused alveolus. This means that the exhaled O2 and CQa tensions must be 
corrected to true or "ideal" alveolar values in order to account for the dilutionary effects of the 
unmetabolized gasses in the deadspace. The "ideal" alveolar partial pressure is atenn used 
herein to indicate that any effects of physiologic deadspace have been incorporated. Physiologic 
deadspace is defined as a region of the lung that is being venfflated but is not being perfused. 

30 The Bohr equation allows correction of the exhaled tensions (of gas X, that is, COj and Oa) for 
physiologic deadspace, if deadspace and expiratory volume; V E can be quantified. 
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V PHY sn/V E = [Pax-Pex]/[Pax-Pk] (17) 

In accordance with the present invention, the physiologic dead space volume, Vphysd, is 
quantified using Fletcher's analysis of the single breath CCfe concentration (Le^ the expired 
volume waveform), which is collected by the cardiac output monitor (see Fig. 5). The exhaled 
gas tidal volume, V & is determined by integrating the expired flowrate using the property- 
compensated pneumotachometer. The gas analyzer of the present invention measures the partial 
pressures in the exhaled gas, Pex, and the inhaled gas, P K . This then leaves the "ideal" alveolar 
partial pressure, Pax, as an unknown. 

Fletcher considered physiological deadspace (V PH ysd) in terms of two components: 
convective airway and alveolar deadspace. Convective airway deadspace (V D aw) extends from 
the lips to the interface between the inspired and alveolar gas. Alveolar deadspace (Vdalv) 
includes deadspace caused by ventflation-perfusion ratio mismatching within terminal respiratory 
units; ventUation-perfusion ratio mismatching between units; venous admixture or right to left 
shunt; and temporal deadspace. Classical methods for es timating Vphysd are based on collection 
of expired air over a period of a few minutes during steady state with simultaneous sampling of 
arterial blood. The partial pressure of CCh is measured from both samples and Vphysd is 
calculated by using Bohr's equation. Fletcher determines this deadspace with a single breath test 

ofC02(SBT-CQa). 

In accordance with Fletcher's technique, a plot of expired CCh fraction (Feoh) against 
expired volume based on the SBT-COi is generated, as shown in Fig. 10. Phase I defines 
conducting airway deadspace. Phase H defines the transition between conducting airway and 
alveolar gas. Phase III defines the alveolar plateau, representing CO2 rich gas from the alveoli 
Phase II and III together represent the effective tidal volume (V Te a) that contributes to gas 
exchange. Fletcher then divides the single breath into regions X, Y, and Z (fig. 10). Area X 
represents tidal elimination of carbon dioxide (Vroca) with effective part of V T . Area Y 
represents wasted ventilation as a result of Vdalv. 

Area Z represents wasted ventilation as a result ofV DAW . The areas are used to express 
the deadspace fractions: Vdaw /Vt — Z/(X + Y + Z); Vdalv /Vt— Y/(X + Y + Z); Vdaw / Vtalv 
- Y/(X + Y), and Vphysd /V t = Y + Z/(X + Y + Z). The dashed line at me top of me graph 
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represents the fraction of the CO2 of a gas in equilibrium with arterial blood (Facoa). The partial 
pressure of CO a in the exhaled breath is not equal to the partial pressure of COa in the alveoli, 
because of admixture of gas from ventiktion-perfosion ratio mismatched alveoli and gas from 
alveoli with efficient gas exchange. If Face* is not known, the value can be estimated by 
5 extrapolatmgtheairveouttoanexhaledvolumeequival^ 

Capacity. Extrapolation is effective for normal individuals andpatients with airway disease, but 

not patients with pulmonary embolism. 

During mechanical ventilation, gas compression increases Vdalv and equipment dead- 
space alters Vtmw» and Phase II may be difficult to distinguish from Phase HL To extract 
10 information from readily obtained data, the present inventors propose a method closely related 
to Boh* concept in which end-tidal carbon dioxide (F^z) is related to me exhaled carbon 
dioxide (Fecc*). VWareaX) is related to me volume of 

breath in which the whole effective volume has CO* fraction Fecc* equal to area ABCDA- The 
ratio, X/areaABCDA = VWEVteff FeoojJ describes the efficiency of the effective volume in 

15 eliminating CO* 

Fundamentally, one needs to know alveolar concentrations in order to calculate 
deadspace,butoneneedstolmo^ To reserve this 

dilemma, the present inventors propose to develop a "virtual patient" consisting of a 
physiological model of respiratory gas uptake and distribution which will be used to mimic the 
real patient, matching the patient's actual measured gas exchange values. One such model is 
aforementioned GUS model, which simulates respiratory gas uptake and distribution using a 
simplified anatomical and physiological model. The simulation uses patient parameters (e*. 
weight, percent body fat), patient physiology (e.g. O2 consumption, pulmonary shunt, systemic 
shunt, cardiac output, deadspace), ventilation mode (e.g., controlled, spontaneous, minute 
ventilation), and gas delivery (e.g. air, O* N 2 0, anesthetics). Predictions of the concentrations 
and uptakes or losses of the respiratory or indicator gasses in twelve (or any appropriate number) 
compartments can be observed. Again, the compartments include eleven tissue compartments 
and one patient breathing circuit of either an anesthesia machine or mechanical ventilator. 

The GUS model is integrated into the computational scheme and is used to predict real 
time rate of exchange of O* and CO2 between the alveolar gas and blood, as well as the exhaled 
partial pressures during normal and nonuniform ventilation, and iterated until what the model 
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output corresponds to what the patient is producing. A specific example of the iterative process 
is outlined in the flowchart shown in Fig. 1 1 and described below. 

First, an initial guess is made for physiologic deadspace by calculating the value by 
Fletcher's technique using the measured COrVolume waveform and the measured end-tidal COi 
5 tension as an initial guess for the alveolar value. From this estimate and Bohr's equation,- 
alveolar O2 is obtained and the procedure outlined below for computing cardiac output is 
performed. The resultant value of cardiac output is provided as an input to the GUS model, and 
a prediction of alveolar and exhaled partial pressures of Cfe and CQjismade. Ifmepaitial 
pressures agree with the measured values then the cardiac output is accepted. Ifthe values differ, 
10 the alveolar CO2 values computed by GUS are used as the next input and the iteration continues. 
When the simulation has converged to within a predeterrnined difference with the measured 
values, cardiac output is accepted and represents the value that has been corrected for shunt and 
deadspace. 

Given the measurements of instantaneous [P B (t), P E (t)] and end-tidal partial pressures 
15 [PetJ of Oa and CO3 and the "ideal" alveolar partial pressures (from Bohr's equation using the 
first estimate of deadspace from the GUS model) for each incremental volume as a function of 
time during the exhaled breath, an equation for the instantaneous respiratory exchange ratio <R) 
is derived from a mass balance of inhaled and exhaled N 2 (which is not taken up), from which, 
both the mixed-venous and arterial values of CQj concentration are computed. From the 
20 defminon of the instantaneous respiratory exchange rauo, 
quotient), 

Rs [C AC o2(t)VV(t)-C^ < 18 > 

After some considerable manipulation, a relationship between the alveolar concentrations of CO2 
and 64, the general alveolar air equation, which includes the effects of all deadspace is derived, 

25 PWt)=(?--P<7)[RFK>2 +F1C02HI -d -R)F K)2 ]-PA02(t)[R+(l -R)F««M1 -0 -R)F«J. 

(19) 

Measured values of Paco* and P A02 are used to determine R and the Pacoj at R = 032 
defines the true mixed venous value. The value at R = Rss defines the true arterial value. Note 
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that the arterial value of C0 2 does not necessarily coincide with the end-tidal value: a cause for 
inaccuracy when using indirect Fick techniques. 

The O2 consumption and CO2 production (V02 and V'coj) for 1he bream is measured. 
Pulmonary capillary blood flow, Qpcbf, is then determined by the indirect Fick technique wim 
5 the assumptions: 1) instantaneous R reflects gas exchange between pulmonary capillary blood 
and alveolar gas and 2) the CCfe dissociation curve is linear, or at least parallel. Using a value 
of, for example, 4.7ml/Iiter per mmHg Pace* for the slope of the whole-body CO2 dissociation 
curve at rest, the pulmonary capillary blood flow, Qpcbf» is: 

QpcBF = V , O 2(Rss-0.32)/[4.7(Pvo2-Paco2)]. (20) 

10 The slope of the CQ2 dissociation curve can be refined if arterid Cb saturation and hemoglobin 
are known or can be adjusted for the metabolic conditions given the approximate level of CQj. 

The pulmonary capillary blood flow approximates the systemic cardiac output when no 
significant pulmonary shunt is present A pulmonary shunt is blood flow that perfuses 
unventilated regions of the lungs. A pulmonary shunt dilutes pulmonary capillary blood flow. 

15 Using conservation of mass, the systemic cardiac output, Qco, can be expressed as: 

Qoo-Qfobf + Qs < 21 > 
where Qs is the pulmonary shunt blood flow. The mass balance for O2 in the blood is, 

C,Qoo = CvQs + CcQpcbf (22) 

where C is the systemic arterial O2 content, Cv is the mixed venous Qj content, and Cc is the 
20 pulmonary capillary O2 content From these relationships an equation for deterniining the flow 
through the pulmonary shunt can be derived, where O2 content can be expressed as aPo2+ 139 
Hb 4 Sato* a is solubility of Q2 in blood, Hb* is hemoglobin leveU Sat is O2 blood satin^on, and, 
P02 is partial pressures of O2 in arterial, capillary or venous blood. The arterial O2 saturation 
(Sataoa) is measured directly using a pulse oximeter. Since the amount of the dissolved O2 in the 
25 arterial, venous and capillary blood (otPoa) is relatively small compared to the amount combined 
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with hemoglobin, it is ignored. The equation for cardiac output thus becomes, 

Qoo = QpcBF[Satco2-Satvoa]/[Satao2-Satvo2]. (23) 

In order to solve the equation, capillary and mixed venous Q2 saturation are obtained from the 
O2 dissociation curve for hemoglobin. The Margaria equation provides a convenient form: 

5 Satoi=[{m-l}K^cs 4 +KPo2{KP^^ C 24 ) 

where, Sates = fractional oxygen saturation, P02 - O2 partial pressure, m = empirical constant 
~ 124, K= empirical constant dependent on pH and temperature. DeFilippe's data shows mat 
the adjustable constant, K, for 6.8 £ pH £7.8 and 293°K £T £ 31 8*K, is, 

K = [0.796{pH - 7.0}/7.0 + 0.152{pH - 7.0>tf.O + 0.021 1 }]exp[-17.72{T - 293}/T]. (25) 

10 The pH effect is especially important because pH varies along the length of the capillary due to 
COz transport. For normal blood (protein content »7%), the relationship of pH to P002 »: 

pH = -03782 In Pcoa +8.806 <26) 

so that me pH nmy be calculated trom me values of 

capillary blood flow. Oncethe various ^contents have been calculated, (^inay be calculated 
15 since(kBFhasalreadybeendctermin^ Qco is the value reported to GUS if the solutions do 

not match and to the clinician when convergence has been achieved. 

While the cardiac output monitoring system of the present invention has been described 

in conjunction with the use of specific Fick techniques and methodologies, it will be understood 

that the system of the present invention can be used with any technique that is based on the Fide 
20 principle and that determines cardiac output from alveolar gas concentrations, including, but not 

limitedtoany: oxygen-based Fick techiuque^ 

carbon dioxide Fick techniques, Fick techniques involving computer modeling using respiratory 
and soluble insert gas Fick techniques. An overview of such techniques is provided by 
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Capek in 'Tick Techniques," Enc of Med. Devices and Instr, Wester, J. D. (editor), Wiley, NY, 
pp. 1302-13 14 (1988), the disclosure of which is incorporated herein by reference in its entirety. 
The particular Fick technique algorithms used to process the measurements provided by the 
cardiac output monitor of the present invention are implemented in software; thus, any of a 
variety of techniques can be performed with the same apparatus simply by running different 
software on the processor. Moreover, software corresponding to a number of different Fick 
techniques can be stored in a memory of the cardiac output monitor, allowing Ihe clinician to 
select a particular Fick technique that is desired or particularly well-suited to the clinical 
situation, or to determine cardiac outout from more than one Fick technique. 

Having described preferred embodiments of a new and improved method and apparatus 
fornon-invasively determining cardiac output and pulmonary function using respired gas analysis 
techniques and physiological modeling, it is believed that other modifications, variations and 
changes will be suggested to those skilled in the art in view of the teachings set forth herein. It 
is therefore to be understood that all such variations, modifications and changes are believed to 
fall within the scope of the present invention as defined by the appended claims. Although 
specific terms are employed herein, they are used in a generic and descriptive sense only and not 
for purposes of limitation. 
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1. A cardiac output monitoring system for non-invasively detennining cardiac output on 
a breath-by-breath basis from respiratory gasses, comprising: 

a flowmeter coupled to configured to receive respiratory gasses and determine respired 
5 flow rate in real time on a breath-by-breath basis; 

a gas analyzer configured to simultaneously determine individual concentrations of plural 
constituents of the respiratory gasses in real time on a breath-by-breath basis; and 

a processor coupled to said flowmeter and said gas analyzer, said processor deteriniiring 
cardiac output on a breath-by-breath basis from the respired flow rate and individual gas 
10 concentrations. 

2. The cardiac output monitoring system of claim 1 , wherein said processor determines 
cardiac output by applying the Fick principle. 
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3. The cardiac output monitoring system of claim 1, wherein said processor accounts far 
effects of deadspace and shunts in determining cardiac output 

4. The cardiac output monitoring system of claim 1, wherein said processor iteratively 
applies a physiological model to determine cardiac output 

5. The cardiac output monitoring system of claim 1, wherein uptake and release of an 
indicator gas is monitored to determine cardiac output 
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